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ABSTRACT 

We present a comparison of halo properties in cosmological simulations of collisionless cold dark mat- 
ter (CDM) and self-interacting dark matter (SIDM) for a range of dark matter cross sections. We find, 
in agreement with various authors, that CDM yields cuspy halos that are too centrally concentrated 
as compared to observations. Conversely, SIDM simulations using a Monte Carlo iV-body technique 
produce halos with significantly reduced central densities and flatter cores with increasing cross section. 
We introduce a concentration parameter based on enclosed mass that we expect will be straightforward 
to determine observationally, unlike that of Navarro, Frenk & White, and provide predictions for SIDM 
and CDM. SIDM also produces more spherical halos than CDM, providing possibly the strongest ob- 
servational test of SIDM. We discuss our findings in relation to various relevant observations as well as 
SIDM simulations of other groups. Taking proper account of simulation limitations, we find that a dark 
matter cross section per unit mass of udm ~ 10^^'^ — 10^^"' cm^GeV^^ is consistent with all current 
observational constraints. 

Subject headings: Dark matter, galaxies: formation, galaxies: halos, methods; n-body simulations 



1. INTRODUCTION 



The cold dark matter (CDM) family of cosmological 
models provides an excellent description of a wide vari- 
l Oty of oboorvational rooulto, from the oarlior obaorvablo 



epochs detected microwa^^e background fluctuations to 



• The number of dwarf galaxies in the Local Group 
is an order of magnitude fewer than predicted by 
CDM simulations, wit h the discrepancy g r owing to- 
wards sm aller masses (Moore et al. 1999a; Klypin et 
al. 199q ). 



' present-day observations of galaxies and large-scale struc 
, ture. A "concordance model" with roughly one-third mat- 
ter and two-thirds vacuum cnergv, either a cosmological 
constant or quintessence (Caldwell, Dave & Steinhardt 



1998), is consistent with almost all cu rrent observations 



on scales > 1 Mpc ( Bahcall et al. 1999| ). 

Recently, improving observations and numerical tech- 
niques have enabled a comparison of CDM scenarios to 
observations on galactic scales of ^ few kpc. The results 



• Hydrodynamic simulations produce galaxy disks 
that are too small and have too little angular mo- 
mentum, yielding a TuUy-Fishcr relation whose zero- 
point is off by several magn itudes from observations 
([Navarro fc Steinmetz 20001). 



• The robustness of rapidly rotating bars in high sur- 
face brightness spiral galaxies ir nplies lower density 



have not been encouraging. There are a number of dis- 



cores than predicted by CDM (de Battista fc Sell- 



wood 1998) 



tinct observations that may be in conflict with predictions 
of CDM: 

The densit^r profile of galaxies in the inner few kilo - 



• Cluster CL 0024-1-1654 is nearly spherical with a 
large, soft core, while CDM t ypicallv predicts tri- 
axial clusters with cuspy cores ( Tyson, Kochanski fc 



partiprs appparg tn hp miirVi s;Viallmynr tVian prpHir-tpH 



DeH'antonio 1998| , though see Miralda-Escude 2000 



by numerical s imulat ions (Navarro, Frenk & White 
1996, hereafter NFW). For density profiles character- 
ized by p{r) oc r" as r ^ 0, CD M predicts a k, —1.5 
with little scatter (Moore et al. 1999), while current 
Ha obs ervations suggest a ~ —0.5 with significant 
scatter J Swaters, Madore fc Trewhella 200(i|; Dalcan- 



ton fc Bernstein 20001, though see Ivan den Bosch fc^ 
Swaters 200^ . 



for a counterexample). 

Each piece of evidence taken individually is perhaps 
not convincing enough to claim that CDM has failed on 
galactic scales. For instance, until recently there was 
controversv amongst sim ulators regarding inner profiles 
( Kravtsov et al . 1998|), but more caref ul simulation s 



have converged on a consistent prediction ( Klypin 2000 ). 
Observationally, inner galactic profiles are uncer tain due 



• The central density of dark matter halos is observed 
to be pc 



to beam smearing effects in H I observations (Swaters 



0.U2 l\lp:pc ' roug hly independent of halo 



Madore fc Trewhella 2000| ; [van den Bosch et al. 2000| ) 



mass (Firmani et al. 2000t), while CDM predicts 
halos with pc <; 1 A/qpc^"^ at dwarf galaxy masses, 
increasing to larger masses (|Moorc et al. 1999). 

-"^Princeton University Observatory, Princeton, NJ 08544 
^Department of Physics, Princeton University, Princeton, NJ 08544 
^Spitzer Fellow 



though samples of high-resolution Ha observations con- 
tinuc to show shallower profil es than predicted by CDM 
( Dalcanton fc Bernstein 2000| ). The number of observed 
Local Group dwarf galaxies may be reconciled with CDM 
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via plausible scenarios for suppressed g alaxy formation 
(e.g. Bullock, Kravtsov & Weinberg 2000), or else compact 
high-velocity clouds could represent the "mi ssing satel 



lites" that are seen in iV-body simulations ( Blitz et al 
1999^] Hydrodynamic simulations of disk galaxy forma- 
tion are fraught with the usual concerns about the effects 
of feedback, artificial viscosity and resolution, though it 
appears the discrepancies above are due to the under lying 
dark matter distribution ( Navarro fc Stcinmetz 2000 ). 

It becomes more interesting to consider alternatives to 
conventional CDM when one recognizes that all these dis- 
crepancies may be symptomatic of a single cause: Dark 
matter halos in CDM simulations appear to be more cen- 
trally concentrated than observed. Recognizing this, var- 
ious authors have recently forwarded a plethora of al- 
ternative dark matter theories that suppress the central 
concentration of dark matter in galaxy halos. Among 



such theories are that the dark matter is warm ( Bommcr 
Larseji &: Dolgov 2000|; |Colin, Avila-Reese & Valcnzucla 



2000 
2000^ 
nov 2|00' 
hat 



I 



Hannestad fc Scherrer 200C ), repulsive (|Goodmai 



fluid ( Peebles 2000|), fuzz y (Hu, Barkana fc Gruzi 



nox 



c lecaying~7Ccn 200C ), annihilating (Kapling- 
Turner 2000p , and the alternative we inves 



tigate here, self-interacting (SIDM; Spergel fc Steinhardt 
2000) . Interestingly, all theories may be tuned to solve the 
problems mentioned above (at least in analytic approxima- 
tions), all theories may be motivated from particle physics 
considerations, and all theories retain the desirable prop- 
erties of CDM on extragalactic scales (though warm dark 
matter is non-trivially c onstrained by this requirement; see 
Narayanan et al. 2000). 

SIDM is governed by a single free parameter, the cross 
section per un it mass apM of the int eracting dark mat- 
ter particle. 



Spergel fc Steinhardt| suggested ctdm 
cm^GeV"^ 



10~^^ — 10"^'' cm^GeV^^ in order to reduce the central 
concentration of galaxy halos by a sufficient amount to 
alleviate the above problems. Intriguingly, this value is 
close to the cross section of ordinary hadrons, motivat- 
ing some particular particle physics candidates for SIDM 
(Steinhardt et al., in preparation). If ctdm is significantly 
smaller than this range, then the optical depth at galac- 
tic densities is much less than unity, implying that SIDM 
would have a negligible effect on the dark matter distribu- 
tion in halos. 

A qualitative picture of the evolution of an SIDM halo 
is as follows: At early times there is no difference between 
SIDM and CDM since the densities and peculiar velocities 
are sufficiently low that collisions are rare; hence SIDM 
makes identical predictions to CDM regarding cosmic mi- 
crowave background fluctuations and the Lyman alpha for- 
est. As the halo forms and grows via gravitational instabil- 
ity, the central density increases. Eventually, collisions are 
so frequent that dark matter particles scatter out of the 
center as fast as they are accreted, and the density growth 
is halted, forming a core. Such a limit is not present in the 
CDM model, where the central density grows unchecked. 
The SIDM core then begins to extend while retaining con- 
stant central density. Heat transfer from the outer parts 
of the halo raises the temperature in the halo core. If the 
halo is truly isolated, then eventually the core thermalizes 
with the exterior resulting in an isothermal halo with a 
steep density profile. This initiates gravothermal collapse. 



where the direction of heat transfer is reversed and the ex- 
terior begins to cool the halo center. However, in a realistic 
cosmological setting, galaxies constantly accrete material, 
keeping the outer halo hot and heat flowing inwards, thus 
delaying core collapse. The interplay between coUisional 
heat transfer and accretion determines whether a halo will 
undergo core collapse in a Hubble time. 

It is important to appreciate that the transport behav- 
ior does not change monotonically with ctdm- For small 
cross-sections, heat transfer increases with udm since the 
frequency of collisions increases; however, for large cross- 
sections, the conductivity n oc cr^^ — > and no heat 
transfer occurs. Thus, as we discuss in the fluid ap- 
proximation is a poor decription of SIDM in the moderate 
cross-section regime proposed by Spergel and Steinhardt. 
Furthermore, the behavior in the moderate cross-section 
regime cannot be surmised by interpolating between the 
fluid and the non-interacting CDM regimes. A proper 
treatment of the SIDM proposal, which includes the inter- 
play of accretion and heat transfer, its non-monotonic de- 
pendence (7dm 7 and the effects of merging demands numer- 
ical simulations designed to explore the moderate cross- 
section regime. 

In this paper we investigate the statistical properties of 
halos in SIDM and CDM in cosmological A^-body simu- 
lations. Our spatial and mass resolutions are sufficient 
to probe the inner regions (~ l/i~^kpc) of small halos 
(^ fewx lO^M©), while maintaining sufficient volume so as 
to have a significant sample of such halos. We us e a Monte 
Carlo technique similar to Kochanek fc Whit^ (2000) to 
model collisions. The primary difference between our sim- 
ulations and prior investigations (discussed in more detail 
in is that we model a cosmologically-significant random 
volume of the universe with self-interaction cross sections 



in the range favored by Spergel fc Steinhardt, enabling 
us to characterize the statistical properties of halos as we 
vary ctdm- 

§0 describes our initial conditions and simulation tech- 
niques using a Monte Carlo A^-body approach. In §^ we 
compare the structural properties of halos in CDM ver- 
sus SIDM models with several cross sections. In par- 
ticular we examine their central densities, inner profile 
slopes a, the mass dependence of a, concentrations, phase 
space densities, and ellipticities, and where possible com- 
pare to observations. In §^ we use lower resolution sim- 
ulations to test the effects of finite particle numbers in 
our Monte Carlo method. In ^ we examine the sub- 
halo population around the largest halo in our simula- 
tions. In §^ we compare our findings to the simulations 
of other groups who have conducted numerical studies of 
self-interacting dark matter, and examine results from a 
wider range of ctdm- We summarize our results and dis- 
cuss observational constraints in §|^. We find that SIDM 
with (TDM ~ 10~^^ — lO"^"* cm^GeV~^ produces halos that 
are in better agreement than coUisionless CDM for a wide 
variety of observations. 

2. SIMULATING SIDM 
2.1. Code and Cosmology 
We use a modified version of GADGET (ISpringel 



Yoshida fc White 200C), a publically-available TreeSPH 



code for distributed-memory parallel machines. Here we 
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only employ the gravitational A^-body portion. We evolve 
a 4 /i~^Mpc randomly-chosen volume of a ACDM universe 
with n = 0.3, Qa = 0.7, Ho = TOkms^^ Mpc"\ and 
(78 = 0.8, similar to the "concordance model" in agree- 



ment with a wide variety of observations ( Bahcall ct al 
1999 |] We generate in itial conditions using COSMICS (|Ma 
& Bei tschinger 1995 ) at z — 49.7, where our particle dis- 
tribution first becomes nonlinear, and evolve to z = 0. We 
employ 128^ dark matter particles in each run, resulting 
in a dark matter particle mass of rUp = 3.6 x IO^Mq, and 
a spline kernel softening of e = l/i~^kpc (i.e. force is New- 
tonian at 2e). To test resolution effects, we also run a suite 
of simulations with 64'^ particles and e — 2h~^kpc. Their 
initial conditions have an initial density field identical to 
the 128"^ runs, constructed by sampling at alternate grid 
points. 

While our 4 h~^Mpc box is small, well below the non- 
linear scale at z = 0, we are interested here in the be- 
havior on scales of a few kpcs, and it is unlikely that the 
missing large-scale power would have a significant effect 
on the inner portions of halos. In addition, our primarily 
conclusions are based on a comparative study between col- 
lisionless and coUisional dark matter for individual halos, 
so we expect these results to be robust to volume effects. 

2.2. Modeling self-interactions 

We have modified GADGET to include self-interactions 
using a Monte Carlo iV-body technique to proba bilisti 



cally incorpora te collisions, along t he same lines as Burk- 
ert (2000) and Kochanek & White, closer to the latter as 
we use Av from individual particles colliding rather than 
setting Av to be the parti cle's velocity; see the discus- 
sion in Kochanek & White. Each pair of particles with 
positions and velocities (ri,vi) and (r2,V2), separated by 
Sx = |ri — r2|/(2e) and dv = |vi — V2I, interact with a 
probability given by 



P /'gcom(f^*^) 



Sv At 



where At is the timestep. 



Xx = 



47r(2e)3 1 
3mp (TDM ' 



(1) 



(2) 



and 



^ J^W{Sx)W{Sx + Sx')diSx') 

Jscomiox) — N— , (3) 

' W^Sx')d{dx') 

where W is the cubic spline kernel used in GADGET. This 
geometrical factor weights the probability of interaction by 
the product of spline kernel-weighted density distributions 
of the two particles at their given separation. The normal- 
ization N is set by requiring that 



fgcom{6x) 47r Sx"^ d{Sx) = 1, 



(4) 



which ensures that when a particle has interacted with all 
its neighbors within 2e, the resulting probability is equiv- 
alent to 

P = ctdm P At, (5) 



where p is the local dark matter density. 

In our code, the scatterings are performed between indi- 
vidual particles at the time that the acceleration between 
those particles is being computed (i.e. during the "tree- 
walk" ) . In order to ensure that all possible scatterings are 
considered, a tree cell is opened whenever it is within 2e 
of a particle, regardless of the opening criterion. 

If two particles scatter, their velocities are randomly re- 
oriented, keeping the magnitudes of their velocities fixed. 
In practice, a running sum is kept of the change in velocity 
due to the interactions that a given particle undergoes on 
every processor, and at the end of the step the velocity 
change for each particle is summed over all processors and 
added to that particle's velocity. In this way, energy and 
momentum are explicitly conserved, even if the scattered 
particles are on different processors, or a particle under- 
goes more than one scatter in a single timestep (which is 
very rare for the cross sections considered here). 

We consider auM = (colhsionless), 10^^'', 
and 10^^^ cm^GeV"^. We also examined (Tdm = 
10-25 cni^GeV"^ and ctdm = 10~^^ cm^GeV'^ in a 64^ 
simulation, which we will examine in §|. The total 
number of collisions per particle in our simulations are 
1.01 for (Tdm = 10"^"' cm^GeV^^, and 6.05 for (Tdm = 
lO^^'^ cm^GeV^^, with slightly lower numbers (0.9 and 
5.3) for the 64^ runs. Note that a factor of ten increase 
in (Tdm translates only to a factor of six increase in the 
n umb er of collisions, since the lowered central densities 
(§ ^.2| ) partially compensates for the increase in (Tdm- All 
runs were performed on Fluffy, a 32-processor Beowulf- 
class machine at Princeton, with each 128'^ run taking ap- 
proximately one week. 



2.3. The Simulated Halo Sample 
We identify dark matter halos using SKID ^ (Spline Ker 



nel Interpo lative DENMAX; see Katz, Weinberg & Hern 



quist 1996), with a linking length of 2e. We only consider 



halos containing 64 or more particles, to ensure a r oughly 



complete sample of such halos in our simulations (Wein- 
berg et al. 1999). Table 1 lists the number of halos for 
identihed in these simulations. 

A specific resolution issue arises from the finite number 
of particles used to probabalistically model collisions in 
the SIDM simulations: The number of particles in a given 
halo must be high enough to properly Monte Carlo sample 
the distribution. As we will show in §0, halos with > 1000 
particles at z = seem to be accurately represented with 
this technique for the simulations considered here. This is 
quite restrictive, but still permits a significant sample of 
halos (roughly 30 in each 128^ run) with which to compute 
statistics. We also use the full sample of halos to examine 
certain aspects, but we will be cautious about interpreta- 
tions made from halos below this "Monte Carlo resolution 
hmit" . 

3. HALO STRUCTURE 

3.1. Halo Profiles 

We determine halo profiles p{r) by spherical averages 
over radii r = e — > 30e, in 20 equal intervals of logr. 



^ ittp: //www-hpcc . astro . Washington. edu/tools/SKID/ 
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Table 1 
Simulation results. 



Model 


(Tdm( cm^ 


GeV-1) 




128^* 








643t 












^'halo 


<-imcd 


Pc.mcd ' 


<A/.mcd 


A'halo 


<-imcd 


Pc.mcd ' 


''A/.mcd 


CDM 







670 


-i.4y 


1.95 


8.0 


111 


-1.58 


1.91 


6.7 


SIDM 


10- 


24 


647 


-0.93 


0.22 


5.6 


106 


-1.44 


0.36 


4.8 


SIDM 


10- 


23 


566 


-0.37 


0.030 


2.6 


89 


-0.53 


0.027 


1.8 


Observed 








« -0.5 


w 0.02 






» -0.5 


« 0.02 





*128'' median values computed for all halos with > 1000 particles (~ 30 per run). 
^64"^ median values computed for all halos with > 500 particles {'^ 10 per run). 

tCentral dark matter density in M0pc-^; values shown axe extrapolated from lh~^ — » 0.5 kpc using amedii 



A sample of 16 halo profiles from our 128^ simulations 
is shown in Figure 1. Each panel shows a halo pro- 
file for (TDM = (solid), 10^^'* (short dashed), and 
10^^'^ cm^GeV^^ (long dashed). Dotted line segments em- 
anating from the innermost radius point of the ctdm = 
curve indicate slopes a = —0.5, —1 and —1.5 for compar- 
ison. The same corresponding halos are chosen from each 
simulation, allowing a casc-by-case comparison of the ef- 
fect of SIDM. The outer halo profiles (r ^ 10/i-^kpc) are 
virtually identical for each halo, showing that the effects of 
self-interactions are limited to the inner few kpcs of halos, 
and confirming that the same halos are being compared in 
the different simulations. The halos in the leftmost column 
are the four most massive ones in our simulations, while 
the halos in other columns are chosen randomly from a 
descending range of masses. Note that rightmost column 
shows halos with roughly 300 particles each, below our 
Monte Carlo resolution requirement of ^ 1000 particles 
(i.e. ;> 3.6 X 1O^M0), thus the effects of self-interactions 
are not necessarily accurately represented in these cases. 

From Figure 1 it is immediately evident that SIDM pro- 
duces halos that have enlarged central cores and shallower 
inner profiles. CDM halos are almost all cuspy (a ^ — 1 

typically), while most (Tdm = 10~^^ cm^GeV"^ cores are 
close to flat, ctdm = 10"^* cm^GeV"^ leads to profiles 
that are intermediate between these two. In some cases, 
non-cuspy CDM halos are seen, especially at lower masses. 
In these cases the halo may have undergone recent merg- 
ing activity that temporarily lowers the central density, 
which is particularly effective in smaller mass halos. Ad- 
ditionally, recent mergers that have not relaxed make it 
difficult to unambiguously identify the halo center about 
which to compute profiles, typically making profiles ap- 
pear shallower. We make no cut in regards to the merging 
history or "isolatedness" of halos, but we do note that the 
missing large-scale power in our simulations will tend to 
generate fewer mergers, and make the largest objects in 
our simulations appear more isolated. 
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Figure 1: 16 selected halo profiles for coUisionless (solid), 
10-^* cm^GeV-^ (dashed), and IQ-^^ cm^GeV-^ (dot- 
dashed). Corresponding halos from each simulation are pre- 
sented, allowing a direct comparison of the effect of SIDM 

on a halo-by-halo basis. The total halo mass in Mq for the 
o"DM ~ halo is shown in the lower right; the SIDM halo masses 
arc typically within 20%. The columns arc ordered by mass, 
with the four highest mass halos shown in the leftmost col- 
umn. Dotted lines from innermost point show reference slopes 
of a = -0.5,-1.0,-1.5. 

At high masses, the effect of SIDM is very promi- 
nent. The upper leftmost halo is Milky Way-sized 

(6 X 10"'^"'^ Mq), and shows a large core of ^ 15ft^^kpc 
for (7dm = 10"^^ cm^GeV"^ (~ 8/i,-^kpc for (Tdm = 
10" cm GeV"^). The core size and difference in inner 
slope become less prominent to lower masses, though this 
could be due to the increasing effects of unrelaxed halos, as 
well as the Monte Carlo resolution issues discussed earlier. 
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We examine these issues quantitatively in §3.4. 

No evidence is seen for SIDM halo profiles that are 
isothermal, as would be expected if the cross section was 
so large that core collapse would occur on timescales sig- 
nificantly shorter tha n a Hubble time. Thi s supports the 
analytic estimates of [Spergel fc Steinhardt that core col- 
lapse on a Hubble time would not occur until ctdm ^ 
10-22 cm^GeV-^ 

Overall, there is a clear trend on a case-by-case basis 
that SIDM results in a reduced central density and shal- 
lower inner slope of the dark matter halo, with increasing 
(TDM having a greater such effect. 

3.2. Central Densities 

Figure 2 shows the central density of dark matter halos 
Pc taken to be the density at our innermost resolved ra- 
dius e, as a function of halo mass. Here we only consider 
halos with more than 1000 particles, where our Monte 
Carlo technique has sufficient numbers to represent the 
collisional behavior (as we will discuss in i^). The central 
halo density of galaxie s is observed to be ~ 0.02 Mq'pc~^ 
( Firmani et al. 2GGGb ), and is consistent with being inde- 
pendent of halo mass. The observed range of halo densities 
is shown as the hatched region, with a majority of their 
data falling towards the lower end of that region. The ar- 
row in the upper left indicates the increase in pc projecting 
the profile from lh~^kpc in to 500 pc, typical of observa- 
tions of dwarf and low surface brightness galaxy central 
densities, using the slope shown. 

SIDM CTDM = 10-23 cm2GeV-^ halos are in good 
agreement with these observations, while ctdm — 
10-2'* cm2GeV-^ produces inner densities that are a few 
times higher, but still marginally consistent with observa- 
tions. In addition, pc in SIDM models show little trend 
with halo mass, in agreement with observations, because 
the core density is set by collisional physics. 




lO"" 10" 

Mh.,o (Mo) 



Figure 2: pc 



i\fhaio for CDM (left panel) and SIDM (right 



10" 



cm^GeV" 



panel). For SIDM, crosses show (tdm = 
while open squares show ctdm = 10^ cm GeV"\ Only halos 
with 1000 or more particles are shown. The hatched region 



indicate s the range of observed pc compiled by Firmani et al 
(2000a). Arrow in upper left indicates how much each value of 
Pc would increase if measured at 500 pc (instead of Ift-^kpc), 
typical of observations, assuming a profile with the slope shown. 

Conversely, the more massive halos in CDM have central 
densities that are too high by at least an order of magni- 
tude already at lft,-^kpc, and because of their cuspy profile 
the disagreement would be much worse at smaller radii, as 
indicated by the arrow in the upper left. Moreover, CDM 



halos have central densities that increase with mass, in 
conflict with observations. 

As such, SIDM halos appear to agree better with obser- 
vations. Table 1 lists the median central density of halos 
with more than 1000 particles in our various models. This 
shows that our simulations reproduce the observed central 
halo densities for ctdm ~ 10^23 cm2GeV-^. 

3.3. Inner Profile Slopes 

We estimate the inner halo profile slope a as the 
slope between the innermost resolved radii, r ^ 1 —^ 
1.5/i^^kpc. Figure 3 shows a histogram of this slope for 
the coUisionless (solid line), 10^24 cm2GcV^^ (dashed), 
and 10- cm GeV-^ (dot-dashed) cases, for all halos that 
have more than 1000 particles. The number of such halos 
in each simulation is indicated in the legend. The quali- 
tative impression from Figure 1 that SIDM produces shal- 
lower inner profiles is quantified in Figure 3. The median 
values of a are indicated by the arrows from the upper 
X-axis, and are listed in Table 1. 

CDM produces halos that have cuspy cores , with amed ^ 
— 1.5. This is consistent with the work of Moore et al 



among others. 25 of the 28 CDM halos have a < —1, indi 
eating that cuspy cores are a common feature of CDM 
models. Conversely, the inner slopes in SIDM models 
are significantly shallower. For ctdm = 10-2^ cm2GeV^^, 
amod ~ —0.9, while for ctdm = 10-2^ cm2GeV-^, amed ~ 
—0.4, with no halos having a < —1. This latter case has a 
median a close to the value preliminarily sugge sted by Ha 
observations of lo w surface brightness galaxies (Dalcanton 
& Bernstein 2000), though a definitive value awaits a more 
thorough analysis of observational biases. 
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Figure 3: Histogram of inner slopes a for coUisionless 
(solid), 10-^'' cm^GeV-^ (dashed), and 10-^^ cm^GeV-^ (dot- 
dashed). Only halos with M > 3.6 x IO^Mq (> 1000 particles) 
are included; the number of such halos is indicated in paranthe- 
ses in the legend. The median values of a are indicated by the 
arrows from the top edge of the plot, and are listed in Table 1. 

The scatter in a is mostly real. There is some scatter 
due to discreteness effects in measuring the inner slope as 
we do. While we could o btain an inner s lope from fit - 
ting a general halo shape (Hernquist 199C; Klypin 2000| ), 
with 5 free parameters the inner slope would be poorly 
constrained by 20 correlated data points, thus we choose 
our simpler definition. Further scatter arises from recently 
merged halos that temporarily have shallower profiles until 
relaxed. However, neither of these effects is very signifi- 
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cant for the large mass halos plotted in Figure 3. Still, we 
choose to quote the median a rather than the mean, in 
order to quantify "typical" halos in these models and re- 
duce sensitivity to outliers, although the mean is similar. 
We note that a significant sca tter in inner slopes is also 
seen i n the observations (e.g. |de Blok, McGaugh fc vaij 
der H filst 1996| ; [Palcanton fc Bernstein 200q ). 



scribed in § B.E 



Our CDM profiles are, at fa ce value, in better agree- 



ment with the analytic profile of Moore et alj (1999) , with 
an asymptotic slope of a = —1.5, rather than an |NFW 
profile having a{r — > 0) = —1. However, profi le fitting is 
a tricky business (as discu ssed in Klypin 2000). By reduc- 
ing the scale radius in the NFW| profile (i.e. increasing the 
concentration), one can push the transition to a slope of 
a = —1 to a radius smaller than l/i~^kpc where we can- 
not resolve the profiles (the "cu sp-core degeneracy"; see 
van den Bosch & Swaters 200C). Thus we suspec t that 



our CDM profiles can also be adequately fit by an NFW 
profile having a large concentration parameter. As such, 
we do not argue for or against either profile form. Our 
simulations can only predict the slope at r « l/i~^kpc, 
and that is what should be compared to observations. 

SIDM appears to be in better agreement with observa- 
tions of the inner slopes of dark halo profiles than CDM. 
At face value, ctdm = 10^^^ cm^GeV"^ is preferred, but 
given uncertainties in observations and simulation tech- 
niques, Ctdm = lO"^'' cm^GeV"^ is probably also consis- 
tent. A similar value of ctdm also reproduces the observed 
central density of galaxies. Such a coincidence is not ex- 
pected a priori, and may represent a significant success of 
the SIDM scenario. 

3.4. Mass Dependence of Inner Slope 

Figure 4 shows a plot of inner halo slope a vs. halo 
mass A/halo for all halos in our CDM (left panel) and 
SIDM (right panel, ctdm — 10^^"^ cm^GeV~^) simulations. 
The curve shows the running median value of a in bins of 
A(logM) = 0.5. 



100 1000 PIG' 
1 FT 




10'2 



Figure 4: a vs. 



(ctdm = 10"^^ cm^GeV"' 



Mhaio for CDM (left panel) and SIDM 
right panel). All halos with 64 or 
more particles are shown. The line shows a running median 
of the a distribution, with a variance computed from all halos 
within each mass bin. 

The CDM case shows almost no trend with mass, with 
the median slope always around —1.2 —1.5. The largest 
halo does have a slightly shallower slope, consistent with 
the trend seen in very high resoluti on CDM simulations of 
individual halos ( Jing fc Suto 200C ) . The scatter increases 
to low mass due to discreteness and merging effects de- 



The SIDM case shows dramatically different behavior, 
suggesting at face value that smaller halos have steeper in- 
ner profiles. However, caution must be used in interpreting 
this result. First, smaller halos have smaller cores, mean- 
ing that a slope measured at a fixed radius (not scaled 
to the halo core size) will result in a steeper slope. For 
CDM this effect is less significant, since the slope remains 
similar from the outer to the inner halo. Second, the 
small numbers of particles in the low-mass halos makes 
the Monte Carlo technique less effective in modeling col- 
lisions, thereby making SIDM appear more like CDM; we 
investigate this issue further in §0. Thus we make no claim 
regarding a trend of a with Mhaio- 

3.5. Mass Concentration Parameter 

As seen in Figure 1, SIDM appears to have the desired 
effect of reducing the concentration of dark matter halos. 
In this section we quantify this effect using a concentration 
parameter, which we define differently than previous au- 
thors in order to facilitate a more direct comparison with 
observations. 

The cano nical d efinition of a concentration parameter 
is given by [NFW as the ratio between the virial radius 
^200 (taken to be the radius at which the halo density is 
200 X t he cos mic mean) and the scale radius of the halo 
in the [NFW| profile. This concentration parameter, how- 
ever, is difficult to compute unambiguously in the case of 
non-isolated halos, and difficult to compare directly to ob- 
servations that seldom extend out to r2oo- Furthe rmore, 
rs is only defined within the context of the specific NFW| 
model, and profile fits are typically degenerate between 
rs and concentration (Klypin 2000). Colin et al. (2000) 



circumvent some of these issues by defining the concen- 
tration as the ratio of the minimum of r2oo and the halo 
radius to the radius that encloses 20% of the halo mass. 
However, this inner radius is dependent on knowing the 
total halo mass, something which is difficult to determine 
observationally. 

Instead, we choose to define a mass concentration pa- 
rameter Cm, based on enclosed mass rather than radii, and 
restrict the scales in our definition to those where obser- 
vations are available, typically r < 20/i~^kpc. We define 



CM = 27 



M(< rin) 
M(< rout)' 



with 



1 



r^out = 8.5 kpc 



220 km s- 



(6) 



(7) 



where Vdrc is the circular velocity of the halo. The choice 
of Tin is arbitrary; here we base it on the Milky Way, as it is 
convenient and results in observationally accessible scales. 
The scaling with i^circ is that expected for self-similar halos 
following the TuUy-Fisher relation. The normalization fac- 
tor of 27 results in a uniform density distribution having a 
mass concentration of unity. A flat rotation curve between 
and rout implies M(r) oc r, resulting in cj\/ = 9. In our 
simulations, we take Wcirc to be the maximum circular ve- 
locity of the halo as output by SKID. 
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Figure 5: Histogram of mass concentrations cm for collision- 
less (solid), 10"^'' cm^GeV"^ (dashed), and 10"^^ cm^GeV"^ 
(dot-dashed). Only halos with M > 3.6 x lO'^Mg (> 1000 par- 
ticles) are included. The median values of cm are indicated 
by the arrows from the top edge of the plot, and are listed in 
Table 1. 

Figure 5 shows a histogram of cm for all halos with more 
than 1000 particles. As expected, there is a clear trend 
for CDM to have more concentrated halos than SIDM, 
with the amount of concentration decreasing with increas- 
ing (TDM- Note that the difference between SIDM with 
10~^^ cm^GeV"^ and 10~^^ cm^GeV"^ is exaggerated rel- 
ative to the difference between the inner slopes of those 
models (cf. Figure 3) . This is because the concentration is 
increased in 10"^'' cm^GeV"^ relative to 10~^^ cm^GeV^^ 
both due to the increased inner slope, as well as the re- 
duced core radius. We also examined the mass dependence 
of Cm and found no obvious trend, but our range of masses 
is small. 

For comparative purposes, our mass concen tration pa- 
rameter Cm may be analytically re lated t o the NFW con- 
centration parameter cnfw- From NFW, 



Figure 6: Mass c oncentration parameter cm, defined in equa- 



te (x) _ 1 In (1 -f CNFwa;) - c^-pwx/ (1 + cnfw^c) 
X In (1 -I- Cnfw) - Cnfw / (1 -|- Cnfw) ] 



(8) 



where x — r'/r2oo, Vc{x) is the circular velocity at x, and 
is the circular velocity at r2oo ■ Our Wcirc is taken to be 
the maximum halo circular velocity, which may obtained 
by maximizing equation ^; this occurs at Xmax ~ 2/ cnfw 
(though we compute it exactly for the results shown be- 
low). 



Let 



„„„ Jci rc = V c{xma.x)lv2w- In appropriate units, 
r2oo (see |NFW| , equation A2). Thus 



W200 



8.5 „ 



(9) 



implying Xin and Xout = Sxin are solely functions of cnfw 
(note that this arises because we defined rjn oc Ucirc)- Us- 
ing M{< r) PC rv'^ir). 



CM = 9 



v1{Xout) ' 



which is purely a function of cnfw- 
tionship is shown in Figure 6. 



(10) 



The resulting rela- 



tion a vs. NFW 



concentration parameter cnfw- 



Figure 6 shows that cm ~ 8, typical of halos in our 
CDM model, corresponds to cnfw ~ 23. This value is 
in agreement with expectations for dwarf galaxies in a 
ACDM model. Conversely, cm ~ 5.6, which is the me- 
dian value for SIDM with lO"^** cm^GeV"^, corresponds 
to Cnf w ~ 11- Note that the minimum value of cm for an 



NFW halo is 3. Thus SIDM with lO"^'^ cm^GeV"\ with 



Cm, med ~ 2.6, produccs halos that typically cannot be de- 
scribed properly by c nfw- T his is because these halos have 
a > —1 typically, so NFW profiles with a{r ^ 0) = — 1 
are a poor fit. This further illustrates why cnfw is a poor 
way to describe halos in general. 

The largest halo in our simulations has a mass com- 
parable to the Milky Way's, « 6 x IO^A/q. The con- 
centrations of this halo are 6.7, 4.4 and 2.3 in CDM, 
o-DM = lO"^'^ cm^GeV"^ and ctdm = lO"^^ cm^GeV-\ 
respectively, cm of the Milky Way halo is somewhat un- 
certain, because of the uncertainty in the rotation curve 
outside the solar circle (i?©) and the effect of baryonic 
mass within i?Q, but we make a rough estimate here. If 
the rotation curve is flat, cm = 9 as stated before. There 
are suggestions that the rotation curve rises somew hat be- 
yond the solar cir cle (though this is uncertain; see Oiling 



& Merrifield 2000), in which case cm is reduced; let us take 



Cm = 8 as a working estimate. The rotation curve mea- 
sures the total mass, so we must correct for the baryons 
to compare with our simulated cm- If we take the frac- 
tion of baryonic mass to be 50% inside Rq and 20% inside 
3i?0, then cm reduces to 5. In addition, baryons adia- 
batically compress the dark matter as they dissipate, so 
we must correct the Milky Way cm further downwards to 
compare to our dissipationless halos. From the analysis of 
Avila- Reese, Firmani & Hernandez (1998), this reduction 
factor is ^ 1.5 — 2, resulting in the Milky Way halo hav- 
ing Cm ^ 3. Thus after reasonable corrections, the Milky 
Way mass concentration appears to be in better agree- 
ment with SIDM than CDM. The rapid rotation of bars 
also suggests a lower concentration for th e Milky Way- 
sized galaxies t han that predicted by CDM ( de Battista & 
Sellwood 199"! ) 



A more direct comparison with simulations may be ob- 
tained from rotation curves of dark matter-dominated, low 
surface brightness galaxies, where baryonic corrections are 
smaller. We expect that this mass concentration measure 
Cm will be relatively straightforward to compute from such 
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rotation curves (e.g. Dalcanton fc Bernstein 2000 ), so we 
look forward to comparisons, cm has the advantage that it 
is independent of halo fitting parameters, as the enclosed 
mass can be obtained directly from the observed circu- 
lar velocity with modest assumptions. In this sense, it is 
a mor e robust comparison than the inner slope and the 
NFW concentration parameter, which are degenerate and 
sensit ive to scales outside those typically observed ( ]van| 
den B asch fc Swaters 2000 ), and the central density, which 
depends on an uncertain contribution from baryons. 

3.6. Phase Space Densities 
A recently popularized measure of the concentrati on of 



dark matter halos is the central phase space density. Dal. 
canton fc Hogan (2000)| find that observed phase space 



densities Q = p/o"^ scale as Q cx a~'^ — (J~'^, where a 
the velocity dispersion, from dwarf sphero idals up to clus- 
ters of galaxies. Observations compiled by Bellwood (2000)| 
suggest a similar relation, albeit with a large scatter, and 
he uses them to arg ue against any form of coUisionless dark 
matter (though see Madsen 2000). 

Figure 7 shows the phase space density Q of dark matter 
within Tin as a function of cr, for halos with cr > 30 km s~^. 
We calculate a as the velocity dispersion around the group 
center of mass velocity, within ri„ (cf . equation |^) . Open 
circles show CDM halos, crosses indicate SIDM (ctdm = 
10~^^ cm^GeV"^) halos. The dashed and dotted lines 
show Q oc a^^ and Q oc cr^^, r espectively, that bracket the 
observations, reproduced from Dalcanton fc Hogan (2000) 
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Figure 7: Phase space density Q vs. velocity dispersion a for 
CDM (open circles), and SIDM with auM = 10"^^ cm^GeV"^ 
(crosses). Dashed and dotted lines bracket observations, show- 
ing scalings of Q oc a^'^ and Q oc cr~*, respectively. 

SIDM generally shows lower phase space densities than 
CDM. SIDM is in somewhat better agreement with obser- 
vations, falling in the middle of the observed range. Mea- 
suring Q in galaxies is a difficult task, because the stel- 
lar velocity dispersion is not necessarily that of the dark 
matter. Furthermore, in rotationally supported galaxies 
the dubious assumption of an isothermal spherical halo is 
used to relate circular velocity to dispersion. Thus Q is 
perhaps not among the most useful observational discrim- 
inants between CDM and SIDM. 

An interesting remark from Figure 7 is that the scal- 
ing of Q(cr ) is roughly the same in bot h models, roughly 



Q (X a . Dalcanton & Hogan (2000) argue that such a 



is not expected based on simple phase packing arguments. 
This further motivates simulations of SIDM that include 
the cosmological growth of halos via dynamical processes 
of merging and accretion. 

3.7. EUipticities 

SIDM produces halos that are more spherical than 
CDM, because of the isotropic nature of the collisions 
(Spergel fc Steinhardt). This is a generic feature of SIDM 



with any significant cross section, since in the inner por- 
tions of halos where collisions are frequent, the velocity 
ellipsoid is quickly isotropized. Thus the shapes of dark 
matter halos provide an important observational discrim- 
inant between CDM and SIDM. 

We compute axis ratios of our halos using the prescrip- 



tion outlined in [Dubinski 
a tensor 



Carlberg (1991)| . They define 



EXiXj 



with a 



^2 



„2 



(11) 



where a is the elliptical radius, and s < q < 1 are the 
axis ratios, and the sum is over all particles with distances 
~ (xi, X2, X3) from the halo center along the axes of the 
ellipsoid. Then, 



M 



and 



(12) 



where M^x > Myy > Mzz are the eigenvalues of A/, q 
represents the axisymmetry of the halo, while s measures 
the halo flattening. Since a depends on q and s, the calcu- 
lation of Mij must be iterated until convergence, which we 
take to be better than 0.01 in q and s. This scheme weights 
particles roughly equally regardless of distance from cen- 
ter, unlike a moment of inertia tensor which weights the 
outskirts heavily, and thus better represents the elliptic- 
ity of the densit y distribution, as shown in Dubinski & 
Carlberg (1991 
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Figure 8: Axisymmetry q and flattening s vs. r for a 6x 10 Mq 
halo (left panel) and a 7x 10^" Mq halo (right panel). Solid lines 
are CDM, dashed lines are SIDM, ctdm = lO"^"* cm^GeV"^ 



0"DM 

and dot-dashed lines are SIDM, udm 
Lower of two curves is s. 
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scaling results from the dynamical assembly of halos, and ctdm 



Figure 8 shows axis ratios as a function of radius 
in our two most massive halos, having masses « 6 x 
10"Mq (left panel) and w 7 x lO^^M© (right panel). 
Sohd hne is the CDM halo, dashed hue is SIDM with 
10-24 cm2GeV"\ and dot-dashed line is SIDM 
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with (Jdm = 10"^'^ cm^GeV^^. q is the upper of two curves 
for a given model. 

CDM halos are fairly triaxial, while SIDM produces ha- 
los that are much closer to spherical. The effect is depen- 
dent on radius, as in the outer regions SIDM and CDM 
become more similar, since the effect of self-interactions is 



I 

confin ed to the inner parts of halos. Still, even at 30fe kpc 

there are significant differences between SIDM and CDM. 
While not stated, this trend with radius is also evident 



of the potential is likely to be dominated by baryons, 
so a comparison to these simulations is not straightfor- 
ward. Further out, perhaps the most direct observations 
of axisymmetry are those for galaxies with HI rings such 
as IC 2006, which suggests a very axisymmetric halo, 
q Ki 0.93 ± 0.08 at k, 13 kpc ( Franx, van Gorkom fc de 
Zecuw 1994 ). Other observations (see packett 199S ) are 



from Figure 1 of |Yoshida et al. (2000b)| 



Figure 9 shows histograms of axis ratios q and s at 
2/i~^kpc (top panels) and lO/i^^kpc (bottom panels) for 
all halos with more than 1000 particles. The median value 
for each model is indicated by the corresponding tick mark 
on the top axis. The difference between CDM and SIDM 
is mo re pronounced at small radii, where CDM produces 



significantly triaxial halos while SIDM halos remain spher- 



more dependent on observational and theoretical uncer- 
tainties such as viewing angle and potential modeling, but 
persistently suggest q ^ 0.8 at 15 - 20 kpc. Both CDM 
and SIDM halos are consistent with these observations. 
Lensing maps of galaxies and clusters offer the best hope 
for mapping the mass potential in the inner halos, which 
should place strong constraints on SIDM. 

Co nversely, observatio ns of s from polar ring galaxies 
(e.g. ^ackett et al. 1994 ) and X-ray isophotes (e.g. Buotc 



ical. At large radii there is a much milder trend to more 
spherical halos with increasing ctdm- s also shows more 
differences than q. 

Figure 9 shows that while CDM produces halos are typi- 
cally more spherical, there is still significant non-sphericity 
in many SIDM halos. In particular, there is a tail in the 
distributions of both CDM and SIDM to smaller axis ra- 
tios. This may be due to asymmetric infall that temporar- 
ily distorts the shape of the density in some halos, par- 
ticularly smaller ones. This also may just be an artifact 
of finite number of collisions in smaller halos. Note that 
the two largest halos shown in Figure 8 show greater dif- 
ferences at 10ft.~^kpc than suggested by the statistics in 
Figure 9. 



& Canizares 1998) suggest a substantial amount of flatten- 
ing, s K, 0.5 ±0.2 at r 15 kpc in the density distribution. 
Such a flattening, if confirmed, may prove troublesome for 
SIDM. The baryonic component would provide a flattened 
contribution, but is not expected to be significant at those 
radii. It is not immediately evident how these discrep- 
ancies may be resolved, but we note that the problem is 
almost as severe for CDM as SIDM in our simulations. 
It is worth mentioning that our small simulation volume 
results in significantly reduced tidal distortion of large ha- 
los, so our simulations may not accurately represent the 
ellipticities of the outer portions of halos. 
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Figure 9: Histogram of axis ratios q (left panels) and s (right 
panels), at 2/i~^kpc (upper panels) and 10/i~^kpc (lower pan- 
els). Only halos with more than 1000 particles are included. 
Tick marks at the upper axis show median values. 

A comparison with observations of halo shapes is as yet 
inconclusive. In the inner portions of dark halos the shape 



4. THE MONTE CARLO RESOLUTION LIMIT 

Our spatial resolution and mass resolution are well- 
understood. However, another resolution issue arises due 
to the Monte Carlo modeling of self-interactions. A Monte 
Carlo method must be sufficiently well sampled, resulting 
in a separate criterion for the number of particles in a halo 
to be well-represented by our simulation technique. In this 
section, we determine this criterion using our suite of lower 
resolution simulations with 64'^ particles described in 

Since we are most concerned with the inner parts of 
halos, we focus on the inner slope as a function of mass 
as the best measure for examining this Monte Carlo res- 
olution limit. Figure 10 shows a plot similar to right 
panel of Figure 4, except for the 64'^ simulation of the 
CTDM = 10"^^ cm^GeV"^ SIDM model. The top axis 
shows the number of particles in these halos. Here we 
compute a at 2/i~^kpc since that is the spatial resolution 
of our 64^ runs. The dashed line is the running median a 
from the 128"^ simulation, computed at 2/i~^kpc. 

Figure 10 shows that for halos with ^ 1000 particles in 
the 64'^ run, the median value of a is within Ict of that of 
the 128'^ run, though consistently lower. By 300 particles, 
the value of a is significantly lower in the 64^ run. The 
reason it is lower is because with few particles, the Monte 
Carlo procedure results in too few interactions to make 
the profile depart significantly from the coUisionless CDM 
case. Thus to lower masses, SIDM looks increasingly like 
CDM when modeled using this technique. 
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Figure 10: a vs. Mhaio for SIDM with ctdm = 
10~ cm GeV-\ for our 64 run. a liere is measured at 
2/i~^icpc, tiie resolution of our 64"^ runs. Solid line shows a 
running median of the a distribution. Dashed line with error 
bars shows a similar curve from 128"^ run. Deviations between 
two at a level > la occur for halos having somewhere between 
300 and 1000 particles. 

Note that this limit is specific to our simulation pa- 
rameters, redshift, and (Tdm, and is not a general state- 
ment about the Monte Carlo Ai"-body technique. The 
limit becomes higher as ctdm is lowered, since collisions be- 
come less frequent, but even for auM = cm^GeV^^, 
1000-f particle halos also appear convergent. We there- 
fore take our Monte Carlo resolution limit to be ^ 1000 
particles. 

We list the values of amed and CM,med for the 64'^ runs 
in Table 1 for comparison with the 128'^ results, where the 
median values here are computed for all halos with more 
that 500 particles in these smaller runs (roughly 10 in each 
simulation). Note that Ofmod is computed at 2/i~^kpc in- 
stead of l/i~^kpc, partially explaining the steeper slopes 
even at the highest masses. In general, the trends indi- 
cated by the 128^ runs are reproduced at this lower reso- 
lution, suggesting that discrete particle effects do not sig- 
nificantly affect our conclusions. We have also examined 
the 128'^ statistics presented previously using a limit of 
300 particles instead of 1000, and our overall conclusions 
remain the same. 
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Figure 11: Subhalo positions within 500/i~^kpc the largest 
object in our volume, in CDM (left panel) and SIDM with 
o"DM = 10"^'' cm^GeV"^ (right panel). Central cross is the 
position of the large halo. 

Figure 11 shows a projected plot of halos within 
500/i~^kpc of the largest halo in our volume, indicated by 
the central cross. Circle sizes are scaled as log(Mhaio), 
with the smallest circles representing halos with M « 
5 X lO^Af0. Left panel shows CDM, right panel shows 
SIDM with CTDM = 10"^^ cm^GeV"^ The positions of 
subhalos are different due to the accumulated differences 
of chaotic orbits within a highly nonlinear potential well. 
Thus a halo-by-halo comparison for these small halos is 
not possible. A careful examination reveals that SIDM 
has fewer subhalos than the CDM distribution, especially 
the smallest ones. 




5. SUBHALO POPULATION 

Self-interacting dark matter is predicted to significantly 
lower the population of subhalos orbiting around large 
halos, thereby bringing simulation predictions into bet- 
ter agreement with observations of the Local Group dwarf 
population. There are two reasons why SIDM has this ef- 
fect: (1) The lowered central concentration and larger core 
radius makes small halos more susceptible to tidal disrup- 
tion, and (2) Dark matter is ram-pressure stripped out of 
small galaxies as they move through the large central halo. 
In our 128^^ simulations, we have one halo that is roughly 
Milky Way sized, having M w 6.7 x lO^M©. In this sec- 
tion we examine the subhalo population around this large 
halo. 



Figure 12: Number of halos within 500/i ^kpc of the 
largest halo, histogrammed by mass, for collisionless 
(solid), 10"^" cm^GeV"^ (dashed), and 10"^^ cm^GeV"^ (dot- 
dashed). Vertical dotted line indicates our 1000-particle Monte 
Carlo resolution limit. 

Figure 12 quantifies this effect, showing the mass func- 
tion of halos within 500/i~^kpc of our largest halo, for 
our three 128'^ simulations. There is a clear trend that 
SIDM suppresses the subhalo population at the small- 
est masses. For 10*^-^ < M < IO^Mq, CDM has 56 
neighboring halos, ctdm = 10^^** cm^GeV^^ has 40, and 
CTDM = 10-^3 cm^GeV"^ has 29. 

While SIDM does reduce the population of smallest ha- 
los, the reduction is not nearly by the order of magni- 
tude required t o obtain agreement w ith Local Group dwarf 
galaxy counts (Moore et al. 1999a). However, the effect 
of subhalo suppression in these SIDM simulations should 
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be regarded as a lower limit to the true strength of the 
effect. The reason is that almost all these subhalos are 
well below our Monte Carlo resolution limit (dotted line 
in Figure 12), therefore their concentrations are approach- 
ing those in CDM models. Thus tidal disruption of these 
halos is not much stronger in SIDM as in CDM, and is 
increasingly similar to lower masses. Furthermore, ram- 
pressure stripping is reduced in effectiveness for the same 
reason that the Monte Carlo technique is less effective in 
these small halos. Hence the numbers of small SIDM ha- 
los are not as significantly suppressed relative to CDM as 
they should be. 

We conclude that the simulations considered here sug- 
gest a weak trend in reducing the number of subhalos with 
increasing ctdm , but due to resolution effects we can make 
no robust quantitative estimates. What is required is to 
simulate a large halo with incredibly high resolution, hav- 
ing subhalos containing thousands of particles to properly 
model the effects of self-interactions. Such a simulation 
is unfortunately beyond the scope of our current compu- 
tational resources. Alternatively, more sophisticated al- 
gorithms are necessary to model self-interactions in small 
halos moving through large ones, which is an avenue we 
are currently pursuing. 

6. COMPARISON WITH PREVIOUS WORK 

A number of authors have investigated SIDM using iV- 
body simulations. The literature divides into two subsets: 
Those that model interactions in the fluid approximation, 
effectively employing a large cross section, and those that 
model self -interactions in the opt ically thin regime as sug- 
gested by ^pergel fc Steinhardt . In both cases, there is 
disagreement over whether SIDM makes halos less or more 
concentrated than CDM. 



resolve these issues. We simply note that the highly op- 
tically thick limit is not the rel evant scenario to test the 



Moore et al. (2000) and Yoshida et al. (2000a) 



late a galaxy cluster within a cosmological context, using 
a treecode with SPH to model interactions in the fluid 
approximation. Both studies resulted in halos that had 
isothermal profiles and were more centrally concentrated 
than CDM halos. This may be because the large effec- 
tive cross section increases heat transfer efficiency, though 
as mentioned in the Introduction, for a sufficiently large 
cross section one expec ts heat transfer to be diminished. 



Yoshida et al. (2000a) suggested that intermediate cross 



sections would likely yield results that were intermediate 
between steep CDM profiles and steeper isothermal pro- 
files, and thereb y argued against SIDM . However, this is 



contradicted by Yoshida et al. (2000b), as well as the re 



suits presented here, confirming that the intermediate case 
results in halos that have long-lived shallow profiles. 

In contrast, Bryan (private communication) uses a adap- 
tive mesh hydrodynamics code to model self-interactions 
in a cosmological volume, and finds that even in the fluid 
limit SIDM produces sizeable, long-lived cores. It is not 
clear why different hydrodynamic codes give different re- 
sults when they should be operating in the same regime. 
Perhaps the effective cross section is larger in the adaptive 
mesh code due to algorithmic differences, reducing heat 
transfer. Another possibility is that a numerical effect in 
SPH in which cold clumps moving thr ough hot halos have 
their drag signi ficantly overestimated ( Tittley, CouchmarJ 
fc Pcarce 1999| ) makes objects rapidly sink into a dense. 



cross secti on range pr oposed by Bpcrgcl fc Steinhardt 



ijurkcrt (2000) and |Kochanek fc Whitc| (2000) simulated 
isolated halos with SIDM having a cross section closer to 
the range of [Spergel fc Steinhardt . They begin with a 
fully formed cuspy galaxy halo and study the evolution af- 
ter interactions are turned on. They both find that halos 
develop a shallow core for some leng th of time, and then 
undergo core collapse. Burkert and [Kochanek fc White 



disagree on the timescales of core collapse; Burkert finds 
tc ~ 16tdyn in a greement with estim ates from two-body 
relaxation, while Kochanek fc White finds a much shorter 



collapse time scale of tc ^ 2tdvn f or the same dark matter 



cross section. Kochanek fc White explain this difference by 
arguing that Burkert| 's method underestimates collisions of 
slow moving particles. We note that our met hod does not 
suffer from this concern, as it is more like Kochanek fc 
Whitd. 



We suggest some possible reasons why the resu lts of 



Kochanek fc White and (to a lesser extent) Burkert are at 



odds with ours. The first is that they begin with a cuspy 
Hernquist profile. This halo evolves rapidly initially (as 
seen in Figure 3 of Kochanek fc White), resulting in an ar- 
tificially large amount of heat transfer. Second, they simu- 
late an isolated halo, ignoring the accretion of dynamically 
hot material during the formation process that would keep 
the outer halo hot and delay core coUape. We note that 
their dimensionless cross section a dm = -^haiofoM/'^s 
converts to ours by a factor of 2 x 10^^'^ cm^GeV"^ 
for Mi,,io » 2 X lO^H lg, and w 5h~^kpc (cf. Figure 1). 
Kochanek fc White's simulation with a dm = 1 already 



isothermal core. It is beyond the scope of this paper to 



produces halos that maintain cores over many dynamical 
times (cf. top panel of their Figure 2). We suggest that 
a somewhat smaller a dm might be consistent with obser- 
vations as well as their limits on core collapse timescales, 
even without considering the effects of accretion and merg- 

Koctianek &: Wtiitc also point out, as we have, that the 
Monte Carlo A^-body technique requires a large number 
of particles for accurate modeling, and show that 10^ par- 
ticles is sufficient. We note that while most of our halos 
do not have that many particles, our largest (Milky Way 
sized) halo has roughly 2 x 10^ particles, and its properties 
are consistent with those o f smaller halos. 

Yoshida et al. (2000b) have now performed a cosmo- 
logical cluster simulation with a cross sections ctdm ~ 
2 X 10-2^ - 2 X 10"23 cm^GeV"^ using a Monte Carlo 
A^-body method. They find shallower central slopes and 
less concentrated cores with SIDM. Our combined results 
span the range from dwarf galaxies to clusters, and are 
in broad agreement with each other if we focus purely 
on the numerical results. For instance, they find that 
cdm ~ 2 X 10"^'^ cm^GeV"^ produces a cluster having 
a core size of 160/i~^kpc. If we extrapolate our largest 
halo's core size (in our ctdm = 10^^'^ cm^GeV^^ run) to 
cluster scales using the expected rcoro ex VcItc, we predict 
a core size of ~ 150/i~^kpc for their cluster. For ctdm ~ 
2 X 10^^** cm^GeV"^, they get 100/i^^kpc core, while we 
would predict ~ 80/i~^kpc for ctdm = 10"^'* cm^GeV""^, 
again in good agreement. They then proceed to scale down 
from clusters to dwarf galaxies using a simplistic argument 
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based on number of collisions, but their scaled result is con- 
tradicted by our direct simulations. We estimate that this 
may be because they use a CDM value for the scale radius 
and cnfw of dwarfs, and compare them to SIDM values 
for their cluster. We find that ctdm = 10~^^ cm^GeV"^ 



observations of the inner portions of dark matter halos 
where CDM appears to fail. In particular, we find: 



produces NFW scale radii that are double that of CDM 
(cf. Figure 6 and discussion); such a factor would go a long 
way towards alleviating the discre pancy. Taking this into 
account, we find the simulations of Yoshida ct al. (2000b)| 
to be broadly consistent with ours. 
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Figure 13; Halo profile of the largest halo in our 64 simula- 
tions, for a range of ctdm values. Halos are progressively less 
concentrated and have larger cores with increasing ctdm. 

In order to explore the high-dDM limit, we ran 64"^ sim- 
ulations of SIDM with (TDM = 10"^^ - 10-22 cm^GeV"^ 
The most illustrative result is to compare the density pro- 
file of the largest halo in all our 64'^ simulations, as shown 
in Figure 13. As seen in Figure 1, there is a smooth 
trend of increasing core radius with ctdm- SIDM with 
= 10-25 cm2GeV-^ is quite similar to CDM, though 
it may also have a core below our 2/i"^kpc resolution limit. 
Increasing ctdm to 10"22 cm2GeV-^, we continue to see no 
evidence for the development of an isothermal core due to 
accelerated heat transfer. The reason is because the colli- 
sions are so frequent in the outer portion of the halo that 
a dense core cannot develop. Instead, collisions randomize 
the dark matter velocities and prevent a smooth radial in- 
flow required to generate a dense core. As dynamically hot 
material accretes onto the halo, heat keeps flowing inward 
and a large core is maintained. Our results are in better 
agre ement with Bryan as o pposed to Moore et al. (2000)| 



and Yoshida et al. (2000a) 



This also illustrates why sim- 
ulating SIDM beginning with an isolated cuspy Hernquist 
profile may not be appropriate for large ctdm; one should 
at least begin with a halo profile that is self-consistently 
stable for a few dynamic times. 

7. SUMMARY 

We present a set of cosmological self-interacting dark 
matter si mulations having cross -sections in the range fa- 
vored by Bpergel fc Stcinhardt (2000). Our simulations 



include the growth of halos from linear fluctuations in a 
random volume of the universe, with sufflcient volume and 
resolution to obtain a statistical sample of galactic halos 
resolved to 1/i-^kpc. We compare the resulting halos on 
a case-by-case basis to those in a collisionless CDM simu- 
lation having the same initial conditions. 

Overall, SIDM is remarkably successful at reproducing 



1. 



The inner slopes of SIDM with ctdm = 
10-23 cni2GeV-^ typical halos have a ~ —0.4 at 
r ^ l/j-^kpc, with some scatter in a. Our CDM 
halos have a « —1.5, in ag reement with previ- 
ous studies (e.g. [Moore et all 1999). SIDM with 
Ctdm = 10-24 cm2GeV- is intermediate between 
these cases, with median a « —0.9. SIDM is in 
better agreement with a preliminary analysis of Ha 
rotation curves of low surfa ce brightness galaxies 
( Dalcanton fc Bernstein 2000 ). 



SIDM with (TDM — 10-23 cni2GeV-^ produces cen- 
tral densities pc ^ 0.01 A/qpc-'^ at l/j-^kpc, and 
shows no trend with halo mass. SIDM with ctdm = 
10-24 cm2GeV-^ has somewhat higher pc values, 
but remains fairly independent of mass. Conversely, 
Pc in CDM halos is much larger than observed, typi- 
cally ^ 0.1 M0pc-^ at lft,-^kpc, and shows a strong 
trend with halo mass. With their steep proflles, 
CDM halos are in signiflcantly worse agreement at 
smaller radii. SIDM is thus is in better agreernent 



with observations, a s has also been argued by Fir- 
mani et al. (2000a 



3. Simulations with SIDM having (Tdm = 
10-24 cm2GeV-^ are intermediate between CDM 
and SIDM with (Tdm = 10-23 cm2GeV-\ indi- 
cating a smooth increase in the effect of SIDM 
with cross section, a result that extends (using 
lower-resolution simulations) from (Tdm = 10-2^ 
10- cm GeV-^ In particular, the generation 
of singular isothermal halos is not seen in any 
of the massive halos simulated, even for (Tdm = 
10-22 cm2GeV-^. This suggests that the dynami- 
cal process of halo growth in a cosmological setting 
helps keep outer regions of halos hot and prevents 
core collapse in a Hubble time. 

4. We introduce a new mass concentration parameter 
Cm based on a more directly observable quantity, 
the enclosed mass within tens of kpc. This halo con- 
centration is significantly lower in SIDM models as 
compared to CDM, providing an observationally ac- 
cessible discriminant that is not dependent on fitting 
a particular profile form. A rough estimate of cm for 
the Milky Way, with large corrections for baryonic 
effects, favors SIDM over CDM. 

5. The central phase space density is lower in SIDM vs. 
CDM mostly due to the reduction in pc- The veloc- 
ity dispersions in the inner regions are quite similar. 
Both SIDM and CDM are consistent with observa- 



tions shown in Dalcanton fc Hogan (2000) , though 
SIDM is mildly favored. 

6. SIDM produces halos that are more spherical, espe- 
cially in their inner regions, as compared to CDM. 
In principle, this is one of the strongest tests of the 
SIDM paradigm, as near the center any value of (Tdm 
that has a non-negligible effect on the dark matter 
distribution will increase the core sphericity, while 
CDM cores are almost always significantly triaxial. 
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However, baryons are likely to dominate the shapes 
of the inner parts of halos, complicating a direct 
comparison, and in the outer parts the differences 
between SIDM and CDM are less pronounced. 

7. The number of subhalos around our largest (Milky 
Way-sized) halo is somewhat reduced with increas- 
ing tJDM, but due to discreteness effects in our 
Monte Carlo A^-body technique, we cannot put ro- 
bust quantitative estimates on the strength of this 
effect. 

Based on these simulations, our currently favored value 
for g pM is somewhere between 10~^^^ cm^ GeV^^ and 
10- 



core sizes would be in better agreement with observations 
if ctdm oc which would result in the effects of self- 

interactions being diminished in hot cluster environments 
as compared to galaxies. Such a scenario occurs naturally 
if the dark matter-dark matter scattering has low-lying 
resonance or bound state contributions, as is the case for 
ordinary nucleons. 

Another theoretical avenue explored in relation to SIDM 
has been modeling the TuUy-Fisher relation. Hydrody- 
namic simulations show that the simulated TuUy-Fisher 
zero point may be brought into agreement with obser- 
vations only if halos are less centrally concentra te d than 
predicted by CDM ( [Navarro fc Steinmetz 2000| ). |Mo & 



Such a cross section simultaneously Mao (2000) determine that a cross section of ctdm 



reproduces both the observed central density and inner 
slope, as well as being consistent with various observa- 
tions considered here, which is non-trivial. In contrast, for 
instance, warm dark matter has difficulty simultaneously 
reproducing the observed central densities, inner slopes 
and s ubhalo population ( Colin, Avila-Reese fc Valcnzucla| 



20001) 



As stated before, inner halo shapes may provide a 
strong discriminant between CDM and SIDM. On galactic 
scales, they are difficult to observe and confused by bary- 
onic contributions. Conversely, clusters provide a cleaner 
test because they have large cores that are not baryon- 
dominated, and their mass distributions are directly ob- 
servable via lensing. Miralda-Escude (2000)| uses the as- 
phericity of cluster MS 2137-23 to (analytically) argue that 
CTDM < IQ-^^-'^ cm^GeV^S effectively ruling out SIDM 
as a solution to halo concentration problems. On the 
other hand, CL 0024-1-1654 is very spher ical, much more 



so than CDM models g enerally predict ( Tyson, Kochan- 



ski & DeH'antonio 1998). Our simulations cannot directly 
address the shapes of clusters, as we have no cluster-sized 
objects in our volume. However, SIDM shows some range 
of halo shapes due to asymmetric infall and unrelaxed mass 
distributions, so it is unclear whether a single object can 
definitively r ule out SIDM. Support f or this statement is 
provided by Yoshida et al. (2000b), whose cluster has 



enough triaxiality to be consistent with MS 2137-23 even 
for CTDM « 2 X 10-24 cm^GeV"^ (their model SlWb), con- 
tradicting Miralda-Escude's scaling argument. We note 



that halo shapes are unaffe cted by annihi lating ( Kapling- 



hat, Knox fc Turner 2000[ ) or decaying (|Cen 20bo| ) dark 



matter, thus they also provide a discriminant between 
these variants and SIDM. 



The cluster core sizes in the simulations of Yoshida et 

30 - 70/i~^kpc 



al. (2000b) are large r than observed {'^ 
Mira lda-Escude 19951: Tyson, Kochanski 



certainly Tor ctdm ' 

(160/i-^kpc), and probably 
10-24 cm^GeV"^ (lOO/i^^kpc). 



T" 



; Dell'antonic 
cm^GcV-i 
2 



X 



X TU 

even for ctdm 
So it may be that SIDM 
has difficulty matching observations at both dwarf galaxy 
and cluster scales. However, adiabatic contraction of 
baryons during the formation of the cD galaxy will re- 
duce the cluster core radius from iV-body predictions, so 
the discrepancy may not be that large. In any case, SIDM 
with (TDM ~ 10^24 _ 10^23 cm^GeV"^ comes remarkably 
close to matching dwarf galaxies, L* galaxies and clus- 
ters given the ^ 10^ range in mass scales, so we reserve 
judgement pending a more careful comparison with ob- 
servations. Yoshida et al. (2000b)| mention that SIDM 



10^23 cm^GeV"^ would produce a correct Tully-Fisher re- 
lation for Ucirc ~ 100 km/s halos. Thus it is conceivable 
that the cross section preferred from halo structure con- 
straints may also alleviate the Tully-Fisher discrepancies. 

Ustriker [2000) suggested that dark matter interactions 
in the centers of halos naturally produce central black 
holes with a mass scaling A/rh c>c cr ' in agreement with 



observations ( Magorrian et al. 1998 ; Fcrrarese fc Merritt 
2000), and he estimates (Tdm <, 10^24 cm^GeV"^ in order 



to avoid central black holes that are too large. However, 
this estimate is based on an p oc profile all the way 
in to the black hole. If such a dense center never arises 
because collisions inhibit its formation, then the limits on 
CDM are weakened considerably. 

Observations of dark matter halos promise to improve 
significantly in the coming years, particularly constraints 
on halo core shapes from lensing and on inner profiles and 
concentrations of halos from Ha rotation curves. If the 
inner parts of dark matter halos are found to be gener- 
ically triaxial, this would be the high place of sacrifice 
for self-interacting dark matter; conversely, spherical ha- 
los would provide strong support for this scenario. The 
main modeling work yet to be is done is an improved ex- 
amination of the subhalo populations in SIDM, as well as 
simulations of a larger range of mass scales. The TV-body 
Monte Carlo approach has difficulty achieving a large dy- 
namic range due to the stringent Monte Carlo resolution 
limit (i.e. discreteness effects in the probabalistic descrip- 
tion of collisions), therefore a different approach may be 
necessary. 

The SIDM simulations presented here are a first at- 
tempt at examining the effect of self-interacting dark mat- 
ter within the context of a realistic halo formation sce- 
nario. The results are quite encouraging that this simple 
variant of the cold dark matter paradigm will alleviate a 
wide range of difficulties faced by CDM on galactic scales. 
We look forward to further investigations and comparisons 
with observations. 
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